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Abstract 

The results of electronic structure calculations for the one-dimensional magnetic 
chain compound Ca3Co206 are presented. The calculations are based on density 
functional theory and the local density approximation and used the augmented 
spherical wave (ASW) method. Our results allow for deeper understanding of recent 
experimental findings. In particular, alternation of Co 3(i low- and high-spin states 
along the characteristic chains is related to differences in the oxygen coordination at 
the inequivalent cobalt sites. Strong hybridization of the d states with the O 2p states 
lays ground for polarization of the latter and the formation of extended localized 
magnetic moments centered at the high-spin sites. In contrast, strong metal-metal 
overlap along the chains gives rise to intrachain ferromagnetic exchange coupling of 
the extended moments via the d-^^2_^2 orbitals of the low-spin cobalt atoms. 
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1 Introduction 



Low-dimensional physical systems are attracting a lot of attention since long 
due to the occurence of exciting physical properties deviating distinctly from 
those known from three-dimensional systems. This includes both, new types 
of ordering phenomena being due to, e.g., Fermi-surface instabilities and new 
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types of excitations as resulting, e.g., from spin-charge separation. Further- 
more, being embedded in three-dimensional crystals, low- dimensional atomic 
arrangements give rise to different length scales and coupling strengths and 
thus may lead to rich phase diagrams. Finally, the existence of magnetic mo- 
ments in such systems raises questions concerning the relevant exchange mech- 
anisms and their anisotropies. 

In recent years, interest especially in one-dimensional systems has focused 
on a new class of transition metal oxides of the general formula AgABOg 
(A' = Ca, Sr, Ba; A, B = transition metal), which crystallize in the trigonal 
K4CdCl6 structure [1]. The space group is R3c (-Dg^, No. 167). In these com- 
pounds, transition metal-oxygen polyhedra form well separated chains running 
parallel to the trigonal axis; see, e.g., Refs. [2,3] for a representation of the 
crystal structure. Space between the chains is filled with the A' cations. Each 
chain consists of alternating, face-sharing AOg trigonal prisms and BOg octa- 
hedra. Typically, the ratio of interchain to intrachain metal-metal distance is of 
the order of two, this fact explaining the pronounced one- dimensionality. The 
chains themselves are arranged on a triangular lattice. As a consequence, in 
addition to showing the abovementioned unique properties these compounds 
allow to study geometric frustration effects as partial disorder and spin-glass 
like behaviour. 

Prominent members of this class are the cobaltates CasCoBOg with B = Co, 
Rh, Ir, which have gained much interest due to their exciting magnetic prop- 
erties. Here we concentrate on Ca3Co206, which has two inequivalent cobalt 
sites, one in octahedral environment (labelled Col) and the other (Co2) cen- 
tered in the trigonal prisms [2,3]. According to powder neutron diffraction 
data reported by Fjellvag et al. room-temperature lattice constants amount 
to a^ex = 9.0793 A and Ch^x = 10.381 A[2]. Data taken at 10 K by Aasland et 
al. resulted in ahex = 9.060 A and Chex = 10.366 A with only slight changes of 
the atomic positional parameters as compared to 298 K [3] . The latter data 
were used in the calculations. 

While the Co-Co distance within the chains amounts to 2.595 A, next-nearest 
neighbour distances across the chains are 5.31 A [2]. This value is somewhat 
larger than the interchain separation of 5.24 A, since neighbouring chains are 
shifted parallel to the trigonal axis by 1/3 of the intrachain separation. Ac- 
cording to the neutron diffraction data, the octahedral Col-0 distances of 
1.916A are considerably shorter than the Co2-0 distances, which amount to 
2.062 A[2]. From the difference it was expected that Ca3Co206 might show 
charge or spin order at low temperatures [2]. Finally, the short Col-Co2 dis- 
tance within the chains, which is a result of the face-sharing of the polyhedra 
and is close to the value of 2.51 A in metallic cobalt, has been taken as indica- 
tive for strong metal- metal bonding [2]. 



2 



Magnetic susceptibility data revealed the onset of magnetic ordering below 
Tci = 24K [3,4,5]. In addition, a second transition was found at Tc2 = 12K 
[4,5]. From Curie- Weiss behaviour above 80 K an effective magnetic moment 
of 5.7 and a paramagnetic Curie temperature of 28 K pointing to predom- 
inant ferromagnetic ordering was inferred. Neutron diffraction measurements 
by Aasland et al. revealed local magnetic moments of 0.08/iB and 3.00/iB for 
octahedral Col and trigonal prismatic Co2, respectively [3,4]. All moments are 
aligned along the trigonal axis. From these data it was furthermore concluded 
that the coupling is ferromagnetic within the chains and antiferromagnetic 
across the chains [3]. While Tci has been associated with the ferromagnetic 
intrachain coupling, Tc2 would correspond to long-range three-dimensional 
ordering [4,5]. 

Magnetization measurements showed a plateau at ~ 1.31//^ per f.u. for low 
field and a steep increase by a factor of three at about 3.5 T, which was inter- 
preted as a ferri- to ferromagnetic transition [3,4]. However, the final magnetic 
moment of about 4/i5 was difficult to reconcile with the values found for the 
local cobalt moments [3,6]. While Aasland et al. proposed a magnetic structure 
with 1/3 of the chains being antiparallel to the remaining 2/3 [3], Kageyama 
et al. argued that at low field and temperatures between 10 and 24 K the 
system assumes a partially disordered antiferromagnetic (PDA) state [4]. The 
existence of this state in different materials is still a matter of controversal dis- 
cussion. While it was also reported for CasCoRhOe [7], recent magnetization 
and heat capacity data allowed to exclude the PDA state for CasCoIrOg [8]. 
For Ca3Co206 it was shown that pressure leads to increase of Tc2 and, hence, 
stabilization of the ferrimagnetic phase [9,10]. Magnetization data taken be- 
low Tc2 on single crystals by Maignan et al. displayed a sequence of up to 
four different magnetization plateaus pointing to a variety of competing mag- 
netic structures [5]. Moreover, AC-x measurements by Maignan et al. showed 
a remarkably strong frequency shift, which has been interpreted as a measure 
of the frustration between different magnetic configurations [5] . These results 
supported interpretations of Ca3Co206 as a triangular lattice of antiferromag- 
netically coupled Ising spins, each formed from a single chain. However, the 
exact nature of the magnetic ground state of this material is still unclear. 

In this letter we report on electronic structure calculations for Ca3Co206, 
which include the spin-degenerate well as spin-polarized configurations 

assuming both the spin-structure proposed by Aasland et al. and the hypo- 
thetical case of ferromagnetic alignment of the chains. Our calculations reveal 
i) low- and high-spin moments, respectively, at the octahedral and trigonal 
prismatic cobalt sites, ii) a rather large oxygen moment due to polarization 
by the high-spin cobalt sites, which together give rise to the formation of 
extended but still well localized CoOe moments and resolve the discrepancy 
between the magnetization and neutron diffraction data, and iii) coupling of 
these extended moments by ferromagnetic exchange via the 3(i states of the 
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low-spin octahedral cobalt atoms. 



2 Methodology 

The calculations were performed using the scalar-relativistic augmented spher- 
ical wave (ASW) method [11,12]. In order to represent the correct shape of the 
crystal potential in the large voids of the open crystal structure, additional 
augmentation spheres were inserted. Optimal augmentation sphere positions 
as well as radii of all spheres were automatically generated by the sphere ge- 
ometry optimization (SGO) algorithm [13]. Self-consistency was achieved by 
an efficient algorithm for convergence acceleration [14]. Brillouin zone sam- 
pling was done using an increased number of k-points ranging from 28 to 408 
points within the irreducible wedge. 



3 Results and Discussion 

3.1 Spin- Degenerate Calculations 

In a first step the electronic properties of Ca3Co206 were calculated with 
spin-degeneracy enforced. The resulting electronic structure and partial den- 
sities of states (DOS) are displayed in Fig. 1. All bands display rather strong 




Fig. 1. Electronic bands and partial DOS of nonmagnetic Ca3Co206. 

dispersion parallel to the chain axis, i.e. along the line F-A, refiecting the one- 
dimensionality of the compound. However, the finite perpendicular dispersion 
points to albeit weak interchain coupling. 

Five groups of bands are identified. In the energy range from -6.4 to -2.0 eV 36 
bands are observed, which trace back mainly to the O 2p states. In contrast. 
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Fig. 2. Partial Col and Co2 3d DOS of nonmagnetic Ca3Co206. Slight broadening 
is due to the DOS calculation scheme [15]. 

the next three groups of bands, which extend from -1.7 eV to the valence band 
maximum, from 0.6 to 1.0 eV, and from 1.8 to 2.2 eV, comprise twelve and two 
times four bands, respectively, and derive mainly from the Co 3d states. Fi- 
nally, bands starting above 4eV are of Co 4s character. As is obvious from 
the partial DOS, p-d hybridization causes substantial d and p contributions, 
respectively, below and above -2 eV reaching about 40% especially in the up- 
per valence band and the conduction bands. Whereas d contributions from the 
octahedrally coordinated Col show up mainly below -3.8 eV and give rise to 
the peak at -4.8 eV, the Co2 d contributions are equally spread over the whole 
energy range between -6 and -2eV. In passing we mention the d admixture 
to the peaks at 0.8 and 2.0 eV, which is of almost pure Co2 and Col char- 
acter, respectively. From the partial DOS the cobalt atoms can be formally 
assigned a d^ configuration, which, together with the absence of charge dis- 
proportionation, is in perfect agreement with the neutron scattering results 
[3]. 

Crystal field splitting due to the octahedral and trigonal prismatic environ- 
ment of Col and Co2 atoms, respectively, is revealed on closer inspection of 
the partial d DOS as displayed in Fig. 2. Note that we used for both sites 
the same global coordinate system with the z axis parallel to the chain axis. 
It deviates from the natural coordinate system for octahedral coordination 
with the Cartesian axes parallel to the metal-oxygen bonds due to rotation 
of the octahedra centered at the Col sites. As a consequence, both the t2g 
and Cg bands arise as admixtures of all five d states (referred to the global 
coordinate system), a situation well known from other trigonal and hexagonal 
systems as, e.g., V2O3 [16]. Nevertheless, at the Col sites, the octahedral crys- 
tal field leads to nearly perfect splitting into fully occupied and empty Cg 
states. In contrast, the trigonal crystal field at the Co2 sites leads to splitting 
into non-degenerate d^^^-r^ as well as doubly degenerate d^y^x'^-y'^ and dxz,yz 
states. While the latter appear almost exclusively in the peak around 0.75 eV, 
the former three orbitals dominate the upper valence band. 
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Apart from the abovementioned strong p-d hybridization two findings are im- 
portant for the understanding of the magnetic properties. First, as has aheady 
been mentioned, the lower and upper conduction band, respectively, centered 
at 0.8 and 2eV are due exclusively to the Co2 d^z^yz and Col Cg states. This 
is related to the fact that the latter states experience cr-type p-d bonding and, 
being antibonding, are thus pushed to higher energies whereas the j>-(i-bonding 
of the dxz,yz orbitals is less strong. In the present situation this difference lets 
us expect that in a spin-polarized calculation the Co2 sites will carry larger 
magnetic moments. Second, we point to the large width of the ^3^2.^2 bands 
as obvious from both the Col and Co2 partial DOS, which is indicative of the 
strong metal-metal overlap along the chain axis. This has been confirmed by 
the calculated crystal orbital overlap population (COOP) and covalent bond 
energies as well as a more detailed analysis of the electronic wave functions. 
According to the latter the d^.z'^-r'^ states give rise to the dispersionless bands 
along F-M-K-F at -1.5 eV and at Ep and show finite dispersion only parallel 
to F-A. Within these bands the bonding states are more Co2-like, whereas the 
antibonding states have increased Col-like character. 



3.2 Spin- Polarized Calculations assuming Ferromagnetic Order 

In a second step, spin-polarized calculations were performed. Yet, in view of 
the predominant infiuence of the ferromagnetic interchain coupling as deduced 
from the experimental data we started out from the hypothetical situation, 
where all chains are aligned in parallel, before turning to the ferrimagnetic 
spin-structure proposed by Aasland et al. [3] . Insofar our procedure is along the 
same line of reasoning as the calculations by Whangbo et al., who, however, did 
not deal with the observed ferrimagnetic structure and concentrated mainly 
on the origin of the cobalt moments [17]. 

For the assumed ferromagnetic order we find a converged solution with well 
localized magnetic moments of 2.73, 0.35, and 0.14 /i^, respectively, on Co2, 
Col, and oxygen. These values refiect the experimental finding of a high and 
low spin state of the cobalt atoms at the trigonal prismatic and octahedral 
sites, respectively, and are close to those reported by Whangbo et al. as well 
as the neutron diffraction results of 3.00 and 0.08//^ for the cobalt atoms. 
The total moment per unit cell, i.e. per two formula units, amounts to exactly 
S.OfiB- Worth mentioning are the rather high magnetic moments associated 
with the oxygen sites. Summing up to almost 1 per formula unit, they pro- 
vide an easy explanation for the discrepancy between the neutron scattering 
data and the magnetization in high field. 

The Co 3d partial DOS given in Fig. 3 display a rather complex relationship 
between the electronic and magnetic properties. Due to the reduced magnetic 
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Fig. 3. Partial Col and Co2 2>d DOS of assumed ferromagnetic Ca3Co206- 




Fig. 4. Partial Col and Co2 3(i integrated moments of assumed ferromagnetic 
Ca3Co206. 

moment the Col "id states show only small spin splitting of the order of 0.3 eV, 
hence, much smaller than the crystal field splitting. As a consequence, the mag- 
netic moment is carried alone by the t2g states. Within this group of bands 
spin-majority ^3^2.^2 states experience a pronounced localization, whereas the 
spin-minority states have a much larger band width. For the remaining 
states things are reversed. Col d contributions below -3.5 eV display small 
spin-splitting but give only a negligible contribution to the magnetic moment. 
This is related to the fact that d-p overlap within the octahedra is strongest 
for states of eg symmetry, which are completely occupied or empty. The situ- 
ation is different at the trigonal prismatic sites, where the spin-splitting is of 
the same order of magnitude as the crystal field splittings. As a consequence, 
spin-majority states are shifted to the energy region, where the oxygen 2p 
states dominate. At the same time, these states experience considerable de- 
localization. In contrast, spin- minority states form sharp peaks at and above 
the Fermi energy but give only minor contributions to the occupied states. 

The different distributions of the electronic states in the octahedra and trigo- 
nal prisms have interesting consequences for the magnetic moments. This be- 
comes clearer from Fig. 4, which displays energy and orbital-dependent spin 
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magnetizations, i.e. the integrated differences of tlie spin-dependent partial 
DOS. For the Col sites it is easily seen that the small magnetic moment is 
mainly carried by the d^z^^r^ states. It starts to build up in the energy re- 
gion above -1.5 eV and is thus due to orbitals well localized within the Col 
atoms. This is completely different for the trigonal prismatic sites. There the 
magnetic moment is rather isotropic, i.e. it grows out of similar amounts from 
all five d states. According to Fig. 4 the d^z^^r^ orbitals contribute 0.6/iB, 
while the doubly degenerate dxz,yz and dx^-y^^^xy states give LS/i^ and O.SfiB, 
respectively. Furthermore, as is clearly seen in Fig. 4, all these moments start 
to build up already well within the energy region, where the oxygen states are 
dominating, thus reflecting the considerable oxygen polarization. To conclude, 
the strong p-d hybridization within the trigonal prisms, which has a strength 
between the rather nonbonding 7r-type and the directed cr-type overlap of the 
octahedra, leads to the formation of a combined isotropic Co2 d-0 p magnetic 
moment. Note that the oxygen p orbitals participate in the local moments and 
not in the exchange interaction between the prisms. This situation goes un- 
der the name extended moment formation and has been previously observed 
in copper oxides [18,19]. Finally, coupling between these extended localized 
moments is mediated by ferromagnetic exchange interaction via the d^z^-r^ 
orbitals of the low-spin Col atoms. 

3.3 Spin-Polarized Calculations with Ferrimagnetic Order 

In a last step, we turn to the ferrimagnetic structure, where 1/3 of the chains 
have their spins antiparallel to the remaining chains. The results obtained from 
these calculations are very similar to those for the hypothetical ferromagnetic 
structure and, hence, provide an a posteriori justiflcation for the validity of 
the latter. Local magnetic moments of 2.72, 0.34, and 0.14 respectively, 
are found for Co2, Col, and oxygen and the total moment per unit cell again 
amounts to ~ S.OOfiB- According to the calculated magnetic moments as well 
as the partial DOS, which likewise look very similar to those obtained for the 
ferromagnetic situation, the differences between the spin up and down chains 
are rather small. Yet, slightly sharper peaks, i.e. stronger localization of the Co 
3d states, are observed for the minority chains. The general situation obtained 
for the ferrimagnetic structure is thus very similar to that grown out of the 
calculations for ferromagnetically coupled chains. 



4 Conclusion 

In conclusion, electronic structure calculations for the magnetic chain com- 
pound Ca3Co206 have revealed strong influence of the local coordination 
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polyhedra on the electronic and magnetic properties. In particular, shorter 
bond lengths and strong cr-type bonding in the octahedra allow for only small 
spin-splitting at these sites. In contrast, due to smaller crystal field splittings 
in trigonal prismatic surrounding the Co2 atoms are found in a high-spin 
state. Strong d-p hybridization causes rather large oxygen magnetic moments, 
which, together with the cobalt moments at the trigonal prismatic sites, take 
part in the formation of extended but still well localized magnetic moments. 
Moreover, intrachain exchange coupling between these effective moments is 
based on metal-metal overlap of the cobalt d3^2_^2 orbitals along the chains 
and mediated by ferromagnetic exchange interaction via the (13^2 _j.2 orbitals of 
the low-spin cobalt atoms. Recently, mapping these results onto a Heisenberg 
model we were able to underline the importance of the interplay of different 
crystal field splittings, d-p hybridizations, and metal-metal overlap for the fer- 
romagnetic intrachain order [20]. Finally, according to a detailed analysis of 
the oxygen magnetic moments as well as the analysis within the Heisenberg 
model the interchain coupling seems to grow out of super-superexchange via 
oxygen states. 
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